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A system of two static quarks, at fixed distances r, and two liglit quarks is studied on an anisotropic lattice. 
Excitations by operators emphasizing quark or gluon degrees of freedom are examined. The maximum entropy 
method is applied in the spectral analysis. These simulations ultimately aim at learning about mechanisms of 
hadronic interaction. 



1. INTRODUCTION 

Features of hadronic interaction can in princi- 
ple be gleaned from two-hadron energy spectra in 
finite boxes 0. For two hadrons containing one 
heavy (static) quark each the relative dis- 

tance r is well defined, and the excitation spec- 
tra as functions of r should give insight into the 
physics of the strong interaction. 

Refined analysis techniques able to deal with 
excitations are desirable. We have employed 
Bayesian inference Q to extract energy spec- 
tra from correlation functions of a set of meson- 
meson operators. In lattice QCD Bayesian tech- 
niques had not been used until only very re- 
cently Experience with this analysis tool is 
thus somewhat limited. In this contribution we 
present selected results from Bayesian curve fit- 
ting with an entropic prior. 

2. OPERATORS AND LATTICE 

The following operators describe systems of two 
pseudo scalar mesons in the 1 = 2 channel 

<^,{t) = J2Sp,s-v (1) 
QAi^thsQAixt) QsiythsqEiyt) 

Up.AA'{xt,yt) U'l.g,g{xt,yt) 

QAixthsgsixt) QB'{yth5qA'iyt) . 

*This material is based upon work supported by the Na- 
tional Scienee Foundation under Grant No. PHY-0073362. 



Heavy and light quark fields are denoted by Q and 
q, respectively, and A, A' , B, B' are color indices. 

Meson fields in $1 are local. Those in $2 are 
spread out over a distance r connected by link 
variable products U p along straight spatial paths 
P between x and y, within time slice t. 

Among the building blocks of the 2x2 corre- 
lation matrix Cij{t,to) = (io)) are light- 
quark propagators, which are computed from ran- 
dom Z2-source estimators living on time slice io 
only, and heavy-quark propagators, which are 
treated in the static approximation. Gaussian 
smearing of the quark fields and APE fuzzing of 
the gauge field links is applied. 

Simulations are done on a L'^ x T = 10'^ x 30 
anisotropic lattice with a (bare) aspect ratio of 
is/ot = 3. We use the tadpole improved gauge 
field action of Q with (3 = 2.4. This corresponds 
to a spatial lattice constant of Os — 0.25fm, aj^ ~ 
SOOMeV. The (quenched) anisotropic Wilson 
fermion action is augmented with a clover term 
limited to spatial planes. Only spatial direc- 
tions are improved with renormalization factors 
Us = (n)^/^, while W( = 1 in the time direction. 
The Wilson hopping parameter k = 0.0679 leads 
to the mass ratio m-^lrfip ~ 0.75. 

3. SPECTRAL DENSITY 

Consider an operator combination = + 
V2^2 and the corresponding time correlation 
function Cy{t,to) = {^l{t)^^{to)). We will fit 
an cj-discretization of the spectral model 

FiPT\t,to) = / dwpT(w)expy(w,i - to) , (3) 
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where expr(a;,i) = e-'^*+"™(-"). The spectral 
density pt{'^) is a discrete sum of (5-peaks 



(4) 



here c„ = e~""^®(~""). The x^-distance of the 
spectral model from the lattice data, involving 
Ct,(ti, io) — ^(ptI^i, io)j is computed with the full 
covariance matrix between all time slices ti ^ Iq. 

We use the discretization Aoj = 0.04, uJk — 
Awfc, k = —50 . . . + 125, and pk = Aw priujk)- 

Assuming minimal information about p the 
Bayesian prior probability is oc e~"^ , where a 
is a parameter and 

•5 = y^(pk - mk - Pk In — ) (5) 

k 

is the entropy relative to a default model m Q]. 
Minimizing the functional 



W[p]^xV^-c^S 



(6) 



then yields the most likely spectral density dis- 
tribution p. Unlike in our implementation is 
based on 



Zw = j [dp] e" 



W 



(7) 



0. 



using simulated annealing, or cooling, /?, 
Local Metropolis updates pk — *■ xpk where a; is 
randomly drawn from the p.d.f. P2{x) = xe~^ 
are employed. 

4. RESULTS 

Analyzing the correlation function of a single 
local heavy-light meson operator, see Fig. re- 
veals two distinct spectral peaks for uj > 0. The 
dominant narrow peak corresponds to the wide 
log-linear stretch of the correlation function. The 
smaller and wider peak originates with data on a 
few early time slices, as closer examination shows. 
This attests to the sensitivity of the method. The 
nature of this excited state is not clear. 

The examples of spectral densities shown in 
Fig. |2 correspond to the two-meson operators (|l|) 
and (|), respectively, both at relative distance 
r = 4. The nonlocal operator data are noisy to an 
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Figure 1. Single- meson time correlation func- 
tion and spectral density with a constant de- 
fault model m = 10~^^ and entropy parameter 
a = 1000. 



extent that precludes the use of standard plateau 
methods. Bayesian inference still yields a spec- 
tral peak, though broad. Empirically, the length 
of a log-linear stretch (plateau) and the width of 
the corresponding peak appear to be inversely re- 
lated. The spectral densities in Fig. || stem from 
averaging over four annealing starts. The spikes 
tend to smooth out for a large number of starts. 

There are three main issues with the current 
approach: Dependence of the results on the en- 
tropy weight parameter a, on the default model 
TO, and on the start configuration for annealing. 

The minimum of the functional W[p] is unique. 
Different annealing start configurations p will 
however move only into the vicinity of the mini- 
mum. (Information about the shallowness of W 
at the minimum may thus be inferred.) It turns 
out that integrated quantities like the peak vol- 
ume and the peak energy 

^„ = ^ Pk ^ \{n\^vitom\^ (8) 
fceA„ 

En = ^n^ X! Pk^k, (9) 

fceA„ 

where A„ is the domain of peak n, are insensitive 
to annealing starts, whereas the micro structure 
(fringes) of the peaks can vary considerably. In 
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this context isolated spikes, like in Fig. ||, are com- 
mon with fine discretization, but have almost no 
effect on and En ■ Those quantities are also ex- 
tremely stable with respect to changes of m and 
a. Typical ranges of stability are a remarkable 
m ~ 10-12 . . . iqO and a ~ 10"^ . . . 10+^ 




Z =165(5) 
E =1.00(3) 
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Zj= 100(11) 
F -1.67(9) 
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Figure 3. Comparison of the $2 versus $1 rela- 
tive excitation strengths \vl{ZiIZ\) and spectral 
energies Ex^i for the ground and excited states as 
functions of the relative distance r. 



creased. Those two mechanisms of interaction 
are repulsive (in / = 2), but move in opposite di- 
rections with r. They possibly compete at small 
r < 1 in the chiral limit j7|. The adiabatic po- 
tential Vadl?') = Min(i?i, i?2) is repulsive for the 
present set of lattice parameters. 

5. CONCLUSION 

Spectroscopic analysis by way of Bayesian in- 
ference is a very powerful method to treat exci- 
tations of hadronic systems. The excitations of a 
heavy-light meson-meson give insight into aspects 
of strong interaction physics. 



Figure 2. Two- meson correlators and spectral 
densities at r = 4. The top and bottom pair 
of figures correspond to $1 and $2, respectively. 



We have used 628 gauge configurations. The 
corresponding statistical errors on are 
comparable to fluctuations due to different an- 
nealing starts. Errors bars in Fig. ^ reflect an- 
nealing start variances. 

Preliminary results for the r dependence of se- 
lected quantities are shown in Fig. ^ The Z2/Z1 
ratio points at an enhancement of $2 -excitations 
(gluon degrees of freedom) over $i-excitations 
(quark d.o.f.) as the relative distance is de- 
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